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The evolution of the Lya forest effective optical depth 
following He ii reionisation 
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ABSTRACT 

Three independent observational studies have now detected a narrow (Az ~ 
0.5) dip centred at z = 3.2 in the otherwise smooth redshift evolution of the 
Lya forest effective optical depth. This feature has previously been interpreted 
as an indirect signature of rapid photo-heating in the IGM during the epoch of 
He II reionisation. We examine this interpretation using a semi-analytic model of inho- 
mogeneous He ii reionisation and high resolution hydrodynamical simulations of the 
Lya forest. We instead find that a rapid (Az ~ 0.2) boost to the IGM temperature 
(AT ~ 10^ K) beginning aX z — 3.4 produces a well understood and generic evolution 
in the Lya effective optical depth, where a sudden reduction in the opacity is fol- 
lowed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low 
density IGM. This behaviour is inconsistent with the narrow feature in the observa- 
tional data. If photo-heating during He ii reionisation is instead extended over several 
redshift units, as recent theoretical studies suggest, then the Lya opacity will evolve 
smoothly with redshift. We conclude that the sharp dip observed in the Lya forest 
effective optical depth is instead most likely due to a narrow peak in the hydrogen 
photo-ionisation rate around z = 3.2, and suggest that it may arise from the modula- 
tion of either reprocessed radiation during He ii reionisation, or the opacity of Lyman 
limit systems. 
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1 INTRODUCTION 



Observations of the He ii iGunn fc PetersonI l|l965l ) troush 
in the spectra of intermediate redshift quasars, coupled 
with large fluctuations in the mean transmission of the 
He n Lya forest, provide the most direct evidence for 
th e tail-end of He ii re i onisation occurring around z = 
3 (IJakobsen et al.l [l994l: iDavidsen et all 1 19961: iHeap et al 



20001 : IZheng et al.H20ol : IShuU et aLri2004l : iFechner et al 
20061 ). However, to date only a handful of ultraviolet 
(UV) quasar spectra suitable for detailed analyses of 
the He ii Lya fore st have been obtained (although see 
ISvphers et al.|[2"009l ). Moreover, the Lya transmission is sen- 
sitive to small He n fractions only, and is thus unsuitable for 
probing the earliest stages of He n reionisation. 

Fortunately, the impact of He ii reionisation on the in- 
tergalactic medium (IGM) can still be probed, albeit in- 
directly, over a wider redshift baseline with the existing 
wealth of high quality, high signal-to-noise H i Lya forest 



data. Firstly, the increased temp erature of the IGM ex- 
pected following H e ii re i onisat ion ("Miralda-Escude fc ReesI 
1994; Abel fc Ha ehneltl Il999l : [Paschos ct al. 2007i) wiU 
ther mally broaden absorpt ion lines in the Lya forest 
( Haehnelt fc Steinmetzjliggsl ). There is some evidence that 
observed line widths are consistent with a sudden increase 
in the IGM temperature around z = 3.3 l|Schave et al.|[2000l : 
iRicotti et~al. 2000), although the error bars on the mea- 
surements are large an d not all studies agree on this result 
([McDonald et al1l200ll ). Secondly, the residual H i fraction 
in the IGM, which is in photo-ionisation equilibrium with 
the metagalactic UV background, will be lowered through 
the temperature dependence of the H ii recombination co- 
efficient (riHi tx T"*^'^). If the temperature of the IGM rises 
suddenly following He ii reionisation, a similarly sudden de- 
crease in the H i fr action, and hence the observed Lya forest 
opacity, w ill result |T hcuns e t al ] |2002l : lBernardi et al ] |2003l : 
iFaucher-Giguere et al.,i2008b '). 
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iBernardi etU] l|2003l ). in a study using 1061 moderate 
resolution quasar spectra obtained from the Sloan Digital 
Sky Survey (SDSS), statistically measured such a depar- 
ture from the otherwise smooth, power law evolution of the 
Lya forest opacity. A sudden decrease of around 10 per cent 
was observed in the effective optical depth at z = 3.2, fol- 
lowed by a recovery to its former pow er law evolution by 

(l2008br) have subsequently 
1)20031 ) result by using 86 



2.9. iFaucher-Giguere et al 



reconfirmed the IBernardi et al 



high resolution, high signal-to-noise spectra obtained with 
Keck/HIRES/ESI and MageU an/MIKE to directly me asure 
the Lya opacity. Furthermore, iDaU'Agho et all ()2008h have 
also recently detected a narrow dip in the Lya forest effective 
optical depth at z — 3.2, albeit at a low level of statistical 
significance (2.6cr), using another set of 40 high resolution 
spe ctra obtained w i th VL T /U VES . 

iTheuns et all l|2002l) (hereafter T02) compared de- 
tailed hydrodynamica l simulations of the Lya forest 
to the IBernardi et al.l l)2003h data and interpreted this 
feature as indirect evidence for photo-heating during 
He 11 reionisation. However, recent theoretical studies in- 
dicate it is difficult to explain this narrow feature in the 
Lya opacity evolution by invoking a rapid temperature 
boost in the IGM alone. (Bolton et al.l (|2009l ) used analyt- 
ical and numerical arguments to demonstrate that there 
are too few hard photons available to heat the entire IGM 
by the required amount (~ 10* K) over the timescale on 
which the observed opacity decreases {Az = 0.1 — 0.2). 
State-of-the-art radiative transfer simulations also indicate 
that He 11 reionisation and reheating will be an extended 
rather than sudden process, with the volume averaged 
IGM temperature ris ing gradually from higher redshifts 
l|McQuinn et al.ll2009l) . A rapid injection of energy into the 
IGM during He 11 reionisation is thus unlikely to explain the 
sharp initiation of the observed dip in the Lya forest opacity. 

Somewhat separate to this argument, however, is the 
explanation for the subsequent recovery of the observed 
Lya opacity to its former power law evolution by z = 2.9. 
Let us assume that the narrow dip in the Lya forest ef- 
fective optical depth is indeed initiated by a large, rapid 
temper ature boost in the IGM following He 11 r eionisation 
(but see lBohon et al.ll2009l : fMcQuinn et al.ll2009l ). Ahhough 
this will not occur over the whole IGM, such rapid heating is 
not necessarily excluded in localised patches close to quasars 
with very hard spectra, where the requisite hard photons are 
abundant. Following reionisation, the thermal evolution of 
the IGM at low density is dominated by the balance be- 
tween phot o-heating and adiabatic cooling due to Hubble 
expansion (|Hui fc Gnedinlll997l : [Theuns et aLlll998l ). How- 
ever, if the Lya opacity evolution were driven by the ther- 
mal state of the IGM alone through the H 11 recombination 
rate, any subsequent recovery in the opacity due to adia- 
batic cooling would occur over a H ubble time - too loni 



to explain the rapid recovery in the iBernardT et all (|200: 
feature. The study of T02 instead demonstrated that this 
interpretation is too simplistic. They found that the hydro- 
dynamical response of the IGM following a sudden reheat- 
ing, coupled with the impact of the extra electrons liberated 
during He 11 reionisation, are enough to drive the Lya opac- 
ity back to its former power law evolution by z = 2.9 as 
observed bv IBernardi et al] (|2003l ) - much earlier than one 
would expect using the simple argument above. 



To date, the only dedicated study of the impact of 
He II reionisation on the Lya forest effective optical depth 
using high resolution hydrodynam ical simulations of the 
IGM has been performed by T02. iMcQuinn et al.1 (|2009l ) 
recently presented a detailed analysis of the evolution of 
the Lya forest opacity using radiative transfer simulations 
of He II reionisation. However, their study necessarily used 
low resolution, post-processed dark matter simulations, and 
they were thus unable to model the impact of the gas hy- 
drodynamics on the recovery of the Lya forest opacity. In 
light of these recent observational and theoretical results we 
therefore re-examine the impact of He 11 reionisation on the 
evolution of the Lya forest opacity using semi-analytic mod- 
elling and high resolution hydrodynamical simulations of the 
IGM. In particular, we shall focus on the explanation for the 
recovery of the dip observed in the Lya forest opacity. This 
paper is therefore closely related to the work presented in 
iBolton et all (|2009l '). where the issues surrounding the initi - 
ation of the feature first detected bv iBernardi et aTl (|2003l ) 
{i.e. the plausibility, or lack of, for the rapid photo-heating 
of the IGM) were examined in detail. 

The structure of this paper is as follows. We begin in §2 
with a brief review of the relationship between the Lya forest 
opacity and the underlying physical properties of the IGM. 
In §3 we use a semi-analytic model to examine the evolu- 
tion of the Lya forest opacity following photo-heating in the 
IGM during inhomogeneous He 11 reionisation. Motivated by 
these results, we then proceed to model the Lya forest opac- 
ity in more detail using hydrodynamical simulations of the 
IGM. The simulations are described in §4, and the evolution 
of the Lya forest opacity in the simulations is presented in 
§5. In §6 we investigate the impact of sudden reheating on 
the IGM gas distribution and peculiar velocity field; we find 
that our simulations are unable to reproduce the narrow fea- 
ture observed in the Lya forest opacity evolution. Finally, 
we consider alternative explanations for the observational 
data in §7 before concluding in §8. 



2 THE FLUCTUATING GUNN-PETERSON 
APPROXIMATION 

The forest of Lya absorption lines observed in the spec- 
tra of high redshift quasars originates from the neu- 
tral hydrogen remaining in th e intervening, low density 
IGM f ollowing H i reionisatio n (iBi et al. 19921: Zhang et al] 
19951: iHernauist et al.l 1 19961 : iMiralda-Escude et al.ril996l : 



Theuns et al 



19981 ). The mean normalised flux of the 
Lya forest, {F) = {/observed //emitted), is the simplest ob- 
servable quantity, and it is often expressed as an effective 
optical depth 



Teff 



-ln(F) = -ln{e-^), 



(1) 



where r is the underlying Lya optical depth in each 
pixel of the spectrum or spectra from which (F) is mea- 
sured. Assuming the IGM is highly ionised and in photo- 
ionisation equilibrium with the metagalactic UV back- 
ground, and the low density IGM (A = p/(p) < 10) fol- 
lows a power-law temperature density rela,t ion, T = TqA^"^ 
iHui fc Gnedinlll997l : Ivalageas et~al]|2002l), the Lya optica l 



depth at z>2 may be written as ( e .g. 
iMcDonald "fc Miralda-Escudell200lh 



Weinberg et~al]|l999l : 
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1.0 



(l+XHc) 



To 



10* K 



0.024 



0.135 



-1/2 



9/2 



A 



2-0.7(7-1) 



(2) 



where fib and ilm are the present day baryon and mat- 
ter densities as a fraction of the critical density, h = 
i?o/100 km s~^ Mpc~^ for the present day Hubble con- 
stant Ho, A = p/{p) is the normalised gas density, Tq is the 
gas temperature at mean density, 7 is the slope of the tem- 
perature density relation and r-12 = Tm/lO"^^ s-i is the 
hydrogen photo-ionisation rate. The power-law temperature 
dependence is due to the case-A H ii recombination coeffi- 
cient, and xho accounts for the extra electrons liberated dur- 
ing He II reionisation; xho ~ 1-08 prior to He ii reionisation 
and xho = 1 .16 a fterwards for a helium fraction by mass 
of F = 0.24 (|01ive fc Sk illman 2004 ). The effective optical 
depth can then be estimated by integrating over all possible 
IGM densities A. 

Eq. ([2J is the Fluctuating Gunn-Peterson Approxima- 
tion (FGFA), and although it ignores the effect of redshift 
space distortions on the Lya forest opacity, it clearly eluci- 
dates the relationship between the opacity and the underly- 
ing physical properties of the IG M. A sudden decrease i n the 
Lya effective optical depth (e.g. iBernardi et al.ll2003l ) can 
thus be attributed to an increase in the IGM temperature, 
either by raising To or changing] 7, an increase in the photo- 
ionisation rate, or a combination of both. On the other hand, 
an increase in the free electron fraction will raise the opacity 
by reducing the recombination timescale. Both an increase 
in the IGM temperature and an 8 per cent change in the free 
electron fraction will be associated with He 11 reionisation. 
Keeping these points in mind, we now proceed to describe a 
semi-analytical model for inhomogeneous He 11 reionisation 
which utilises this useful approximation. 



3 THE Lya FOREST OPACITY DURING 
INHOMOGENEOUS He 11 REIONISATION 

3.1 A semi-analytic model 

Before analysing our detailed hydrodynamical simulations, 
we first consider a simplified semi-analytic model which illus- 
trates th e difficulty of reproduc ing the localised feature ob- 
served bv lBernardi et alTl|2003l ) in the Lya forest. Although 
this model is unable to compute many of the detailed prop- 
erties of the Lya forest (such as line broadening and the 



^ Lowering (raising) the value of 7 will increase the temperature 
in underdense (overdense) regions of the IGM while decreasing the 
temperature in overdense (underdense) regions. Whether raising 
or lowering 7 subsequently decreases the effective optical depth of 
the Lya forest thus depends on redshift. At 2 > 3 the Lya forest 
opacity is mainly sensitive to underdense regions in the IGM. 
Lowering 7 therefore decreases the Lyo effective optical depth 
at 2: > 3 by producing hotter, more highly ion ised voids (e.g. 
iBolton et al]|2005l : IPaucher-Giguere et al.^l2008a^ . Typical values 
adopted for 7 in Lya forest models lie in the range 1 < 7 < 
1.6, although there is some evidence to suggest a more complex 
relationship between temperature and density, perhaps with 7 < 
1, may be r equired to reproduce the observed Ly a forest flux 
distribution llBecker et ahlboOTl : IBolton et al.ll2008() . 



impact of peculiar velocities), it has the virtue of includ - 
ing inhomogeneous reionisation (|Furlanetto fc Oh! I2OO8IJ ) , 
which our hydrodynamical simulations d o not. The model 
is base d on the calculations presented in IFurlanetto fc OhI 
(|2008al ). who examined how He 11 reionisation can affect the 
temperature-density relation of the IGM, and we refer the 
reader there for more details. 

In brief, the model has three parts. The first determines 
the reionis ation history of gas elem ents of a given density in 
the IGM (IFurlanetto fc OhI l2008al ') . For most of our calcu- 
lations, we assume that the probability that an element is 
reionised at any particular time is independent of its den- 
sity, simply tracing the mean ionisation history, Xi{z). We 
also compare to a "density-driven" model in which high- 
density regions are ionised first (because they sit near the 
biased regions that host quasars). However, recent numeri- 
cal simulations of He 11 reionisation suggest that the ionised 
region s are large and rare enough that this correlation is 
weak (|McQuinn et al.ir2009l ). so we usually use the density- 
independent model. We will assume that He 11 reionisation 
ends at zhc ~ 3.2, consistent with the calculations below, 
and that hydrogen reionisation (which only affects the tem- 
perature of gas for which helium is still singly-ionised) occurs 
at zh = 8; the latter has little effect on our calculation. 

To make this comparison as straightforward as possible, 
we assume that Xi{z) is simply proportional to the mass in 
galaxies with m>5 x 10^^ M© (in other words, these mas- 
sive galaxies contain supermassive black holes that have 
gone through bright quasar phases). This provides a rea- 
sonable approximation to the qua sar emissivity (comp ared 
to, e.g., the luminosity function of lHopkins et al.ll2007l ) and 
also leads to fast He 11 reionisation, with over 70% of the 
ionisation occurring after z = 4.2, i f He 11 reionisation ends 
at z = 3.2 (IFurlanetto fc Ohll2008bl '). 

The second ingredient is to follow the thermal evolu- 
tion of each gas element after reionisation . We ag ain use 
the method present ed in Furlanetto fc ^ (|2008al ). which 
is in turn based on iHui fc GnedinI (|l997D . We include all 
of the relevant atomic cooling, heating, recombination, and 
photo-ionisation processes, as well as adiabatic expansion 
(and collapse for overdense regions). After reionisation, 
the dominant processes are photo-heating and adiabatic 
cooling, which together determine the asy mptotic, nearly 
powe r law, tem perature-density relation (|Hui fc GnedinI 
il997i : iHui fc Haiman. 200311 . 

The chief uncertainty in this model is the initial tem- 
perature after reionisation, T, which depends in a non- 
trivial man ner on the spectrum of the photons ionising the 
gas parcel (|Abel fc Haeh nclt 1999!: 
Paschos et al. I I2OO7I : IBolton et al 



Tittlev fc Mciksin 2007|; 



20091 : .McQuinn et al] 



20091) . This may itself be inhomogeneous in the IGM, be- 
cause low-energy photons will be absorbed near their hosts, 
leaving many of the ion isations to be done by hard photons 
(|Abel fc Haehneltlll999h . We therefore consider a range of 
possibilities for T. Our fiducial model takes T = 4 x 10^ K, 
larger than is typically expected in order to exaggerate the 
effect of heating {e.g., IBolton et al.l I2OO9I : iMcQuinn et al.l 
2009). 

Finally, in order to compute the average optical depth 
we need to assume an IGM density distribution, Pv(A) (av- 
eraged by volume), as well as the distribution of tempera- 
tures and ionisation rates at each density. We will use the 
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densit y distribution recommended bv lMiralda-Escude et alj 
l|2000l ): 



(A 



-2/3 



2(25o/3)2 



dA. 



(3) 



This form fits cosmological simulations at z — 2-4 quite well. 
Note that the underlying simulation had somewhat different 
cosmological parameters than the currently preferred values; 
it is, however, accurate enough for the qualitative calcula- 
tions that follow. 



3.2 The Lya effective optical depth 

To compute the Lya effective optical depth, we assign each 
volume element of known density an optical depth accord- 
ing to Eq. ((2)0. Our model makes no predictions about the 
amplitude of the ionising background, and for simplicity we 
set r_i2 = 1 over the entire redshift range z = 2-5. This 
is consistent with constraints derived by comparing obser- 
vational data to sim ulations of the Lya forest at 2 < 2 < 4 
l|Bolton et aklbOOSl ). and yields a mean transmission reason- 
ably close to the observed values. 

We choose the temperatures by following gas elements 
of the appropriate density after their most recent reionisa- 
tion event using the thermal evolution code described above, 
with the reionisation redshifts distributed according to the 
overall ionisation history. We then calculate 



dAPv{A) / drP(r|A)e' 



r(T,A) 



(4) 



where P(rjA) is the probability distribution of tempera- 
tures for elements at density A and TcS is the effective op- 
tical depth. Again, we emphasise that this ignores peculiar 
velocities of the gas, line blending (and indeed the wings of 
every line), and the clustering of the absorbers, but it pro- 
vides a qualitative description of the evolving transmission. 

Fig. [1] shows some example histories computed in this 
manner. As a basis for comparison, the dotted curve ignores 
Ife 11 reionisation and assumes a constant ionising back- 
ground and IGM temperature (T = 20,000 K). This does 
not quantitatively match the observed evolution; it merely 
serves to show that, without He 11 reionisation, our model 
produces smooth, nearly featureless evolution. 

The other curves assume that helium reionisation ends 
at zne = 3.2. The short-dashed, solid, and dot-dashed curves 
assume Ti = (3, 4, 6) x fO'* K, respectively (all relatively 
large in order to exaggera te the effect; see iBolton et al.l 
I2OO9I : iMcQ uinn et al] |2009| for self-consistent estimates of 
the plausible temperature range). The long-dashed curve 
takes Ti = 4x10"^ K and assumes that reionisation is density- 
dependent (iFurlanetto fc Oh 2008a ): this aspect of reionisa- 
tion has only a small effect on the mean transmission. This 
is because the Lya forest is sensitive to only a narrow range 
of densities near the mean, where the density modulation is 
modest anyway. 

According to these models. He 11 reionisation can cer- 
tainly induce a feature in T^g, so long as the temperature 

^ We drop the assumption of a power law te mperature-density 
relation in this instance; see iFurlanetto fc Oh. 1 2008a ') for some 
example temperature distributions. 




Figure 1. Mean transmission histories in our semi-analytic mod- 
els. The dotted curve shows T^f{{z) if we ignore He II reionisation 
and set T = 20, 000 K throughout. The other curves include 
He II reionisation at zjjc = 3.2, as described in the text. The 
short-dashed, solid, and dot-dashed curves assume Ti = (3, 4, 6) X 
10"* K, respectively. The long-dashed curve assumes T = 4x 10'* K 
but that reionisation proceeds from high to low densities. 



increase is large enough. In all cases, Teff falls relatively 
steeply until zho before levelling off and returning closer to 
the expected evolution without He 11 reionisation. The shape 
is generic within these reionisation models, although note 
that the downward turn in is still signific antly shallower 
than t hat observed in the observation al data l|Bernardi et al.l 
l2003i : iFaucher-Giguere et al.l l2008bl ). The pre-reionisation 
phase can be steepened by making reionisation occur faster, 
but the model assumed here is not far from empirical esti- 
mates of the evolution of the quasar emissivity. Moreover, 
compressing He 11 reionisation into a short time interval lim- 
its the temperature jump that it can induce, because the 
photo-ionisation timescale associated with high-energy pho- 
tons (which provide the most efficient heating) is quite long 
(|Bohon et al.ll2009l ). 

However, the post-reionisation evolution is completely 
generic and easy to understand. In this regime, the thermal 
evolution is dominated by the competition between photo- 
heating and cooling by adiabatic expansion. The timescale 
for the cooling is therefore the expansion, or Hubble, time 
- a substantial fraction of the age of the Universe. There 
is no way to avoid this behaviour for the temperature his- 
tory, and to the extent that the Lya forest depends only on 
these temperatures. He 11 reionisation cannot induce a nar- 
row feature in the transmission which recovers quickly to its 
pre-reionisation evolution. 

Of course, we have emphasised that the forest is indeed 
more complex than this model, because of peculiar veloci- 
ties, geometric effects, and line blending. In principle, these 



ap pe 

an effect to explain the iBernardi et al.l l|2003l ) feature. In 
their simulations, once the reheated gas in the IGM became 
overpressurised with respect to its surroundings it started to 
expand, resulting in a sudden change in the peculiar veloc- 
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ity gradients in the IGM. This extended the Lya absorption 
lines in redshift space by shifting absorption from the satu- 
rated hne cores to the wings, increasing the hne equivalent 
widths and hence the mean Lya opacity. In the remainder 
of this paper, we will examine all of these effects in more de- 
tail using new hydrodynamical simulations that do include 
a detailed reconstruction of the Lya forest. 



4 HYDRODYNAMICAL SIMULATIONS OF 
THE Lya FOREST 

4.1 Initial conditions 

The hydrodynamical simulations in this study are performed 
using an upgraded version o f the publicly available parallel 
Tree-SPH code GADGET-2 (|Springelll2005h . All simulations 
have a box size 15h~^ comoving Mpc and contain 2 x 400^ 
gas and dark matter particles. The mass of each gas parti- 
cle is 9.4 X 10^ M0 and the gravitational softening length 
IS 1/30* of the mean linear interparticle spac ing. This ade- 
quately resolves the Lya forest at 2 < 2 < 4 (|Theuns et al.l 
1 19981 : iBolton et al.l [20081 ) and provides just over twice the 
mass resolution of the simulations used in the T02 studyQ 
Star formation is included using a simplified prescription 
which converts all gas particles with overdensity A > 10^ 
and temperature T < 10^ K into coUisionless stars, signifi- 
cantly speeding up the simulations. Outputs are saved every 
Az = 0.05 in the redshift range 2 < z < 4, enabling a very 
fine sampling of the simulation data with redshift. 

The ionisation state of the gas particles is computed 
in the optically thin limit using a non-equilibrium ionisa- 
tion algorithm which follows the abund ances of six species: 
H I, H II, He I, He 11, He iii and e~ (|Bolton fc Haehneltl 



200781). The ult r aviole t background (UVB) model of 



Haardt fc Madaul ()200lh (hereafter HMOl) including emis- 



sion from quasars and galaxies is used to compute the 
H I, He I and He 11 photo-ionisation and heating rates. This 
model is in good agreement with observationa l constraints 
on the metagalactic H i photo-ionisation rate (|Tvtler et al.l 
l2004l : lBolton et al.llioosl 'l. ahhough we set the HMOl photo- 
ionisation and heating rates for all species to zero at z > 6 
in order to match the simulations of T02. 

The UVB is assumed to be spatially uniform in all our 
simulations. This will be a reasonable approximation at the 
H I photo-ionisation edge at 2 < 4, when the mean free path 
for hydrogen ionising photons is much l arger than the aver- 
age separation between ionising sources llBolton fc Haehneltl 
l2007bl : iFaucher-Giguere et al.l l2008al : iFurlanettol l2008h . 
However, a full radiative transfer implementation is required 
to model the impact of inhomogcneous He 11 reionisation 
on the Lya forest opacity (e.g. Mas c Ui fc Ferrara 20051: 



Tittley fc MeiksinI [ioOTi: iPaschos et al.l boOTi : iBohon et al.l 
2OO9I : iMcQuinn et all l2009l ). Although He 11 reionisation 



does not directly impact Lya forest effective optical depth, 
our simulations do not capture the effect of the resulting fiuc- 
tuations in the IGM temperature, which do impact on the ef- 
fective optical depth through the temperature dependence of 



^ T02 use a modified version of Hydra l lCouchman et aTlll995l) to 
run fiydrodynamical simulations in a 12h~^ comoving Mpc box 
with a gas particle mass of 2.0 X 10^ Mq. 



the H II recombination coefficient. However, our simulation 
volume is 15h~^ comoving Mpc, which is small on compar- 
ison to the ~ 50 comoving Mpc scal es on which large tem - 
perature fluctuations are expected (|McQ umn et al.l l2009l ). 
Although not ideal, our assumption of instantaneous reheat- 
ing on this smaller scale is therefore a reasonable approxi- 
mation. Furthermore, note that HMOl implicitly assumes 
that He 11 is reionised at all redshifts, i.e. He 11 only exists 
in the dense Lya absorbers which they model. Radiative 
transfer effects obviously modify the for m of the UV back- 
groun d prior to He 11 reionisation (e.g. iMadau fc Haardt 
120091 ) ■ both in He II ionising radiation above 4 Ry and in re- 
processed radiation from He 11 Lyman series and two photon 
emission. Thus, the UV background could be quite different 
before reionisation is complete. We comment on this more 
in ^ 

The simulations were all started at 2 = 99, with 
ini tial conditions g e nerat ed using the transfer function 
of lEisenstein fc hJ (|l999l ). The cosmological parameters 
adopted are = 0.26, = 0.74, Qy^h^ = 0.024, h = 0.72, 
(jg = 0.85 and ris = 0.95. These are consistent with the fifth 
year Wilkinson Micro wave Anisotropy Probe (WMAP) data 
iDunklev et al.ll2009^ ■ aside from a slightly larger value for 
the matter power spectrum normalisation. The gas is as- 
sumed to be of primordia l composition with a he lium mass 
fraction of Y = 0.24 (e.g. lOlive fc SkiUmanlliooil ) . 



4.2 Thermal and ionisation histories 

Five hydrodynamical simulations, listed in Table [T] were 
performed for this study. The simulations use different mod- 
els for the thermal and ionisation evolution of the IGM, but 
aside from the S3 model, which we shall discuss later, are 
identical in all other respects. The different thermal histo- 
ries are constructed by increasing the HMOl He 11 photo- 
heati n g rate in the s imulations {e.g. T02; iBolton et al.l 
I2OO5I : iJena et all |2005| ). This mimics spectral hardening 
due to radiative transfer effects during He 11 reionisation by 
boosting the moan excess en ergy per H e 11 photo-ionisation 
(Abel fc Ha ehnelt 199_9; B olton et al.|[2004 ). Note that we 
only model the jump in heating rates, which affects the ther- 
mal evolution; apart from model S2 which we discuss below, 
we do not model the jump in the He 11 fraction itself. 

The IGM temperatures at mean density. To, in all five 
simulations are displayed in the left panel of Fig. [5] as a 
function of redshift. These dat a are compar e d to obser- 
vational constraints obtained by ISchave et al.l l|2000l ) from 
an analysis of absorption line widths in the Lya forest. 
Note that these thermal histories are not fits to the ob- 
servational data; they are instead merely chosen to be rep- 
resentative of IGM thermal histories presented in the lit- 
erature. The SI model (solid curve) closely resembles the 
To evolution in the simulations of T02, with a sudden in- 
crease {Az — 0.1) in temperature, AT ~ 10* K, beginning 
at z — 3.4. Note that this temperature boost occurs over 
a much shorter timescale th an that in our semi-ana lytical 
model of He 11 reionisation (|Furlanetto fc "Ohll2008a'). The 
El model (dot-dashed curve) is qualitatively similar to re- 
cent results from detailed three dime nsional radiative trans - 
fer simulations of He 11 reionisation (|McQ umn et al.ll2009l ). 
The temperature boost develops over a longer timescale as 
quasars gradually photo-heat the IGM. The third simula- 
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Figure 2. Left: The temperature of the IGM at mean density, To, in the five simulations used in this study: SI (soUd curve), S2 (long-dashed 
curve), S3 (short-dashed curve), El (dot-dashed curve) and Nl (dotted curve). The temperatures are given at i nterva ls of Az = 0.05 over 
the redshift range 2 < 2 < 4. These data are compared to observational constraints obtained bv lSchave et al.l 1I2OOOI ) from an analysis of 
absorption line widths in the Lya forest. Centre: The corresponding volume weighted He II (upper panel) and H I (lower panel) fractions. 
Right: The slope of the power law temperature-density relation, T = TqA''"^. The values are obtained by a least squares power law fit 
to the volume weighted temperat ure-density plan e in the simulations at 0.1 < A < 1 and T < 10^ K. The data are again compared to 
observational constraints from iSchave et al.l (|2000|). Note that the SI and S3 models are almost indistinguishable in all three panels. 



Table 1 . Hydrodynamical simulations used in this study. All the 
simulations have a box size of 15h~^ comoving Mpc and contain 
2 X 400"^ gas and dark matter particles. 



Model Thermal history description 



81 Sharp temperature boost 
El Extended temperature boost 

Nl No temperature boost; control model 

82 Similar to 81, but with a rapid change in rie 

83 Identical to 81, but with a stricter timestep limit 



tion, Nl (dotted curve), is included as a control model, and 
has no temperature boost. 

The corresponding volume weighted He 11 and 
H I fractions in the simulations are displayed in the 
central panel of Fig. [S] In the SI, El and Nl simulations, 
both H I and He 11 reionisation commence at 2: = 6. This 
choice is deliberate; these three models are designed to 
exclude the 8 per cent increase in the free electron fraction 
following He 11 reionisation (see discussion in §2). Instead, 
they will be used to explore the impact of differences in 
the IGM thermal state alone on the Tcff evolution. Sudden 
He II reionisation ai z = 3.4 is instead included in a fourth 
model, S2 (long-dashed curve), which is tailored to have a 
similar temperature at mean density to the SI model. The 
S2 simulation is therefore the most similar to the model 
used by T02, who also assumed H i reionisation ai z — 6 
and He 11 reionisation at z = 3.4. 

The right hand panel in Fig. [2] displays the evolution of 
the slope of the temperature-density relati on, T = TqA^"^ 
jHui fc GnedinI 1 19971 : IValageas et al] |2002| ). in each of the 
simulations. The power-law index 7 gradually increases to- 
wards lower redshifts in the El and SI models, asymptot- 



ically approaching the upper limit of 7 ~ 1.6 achieved by 
the balance between photo-heating and adiabatic cooling in 
the low density IGM. Note however, the SI and S2 mod- 
els, although having very similar To values, exhibit different 
behaviour for 7. This is due to the different ionisation his- 
tories adopted in the two models. During the reionisation of 
He II at 2 = 3.4 in the S2 model, photo-ionisation equilib- 
rium is no longer a good assumption and the He 11 photo - 
heating rate is independent of density (iBolton et al.ll2009l ). 
This is because the He 11 fraction - /hoII = »^HeII/'^He ~ 1 
immediately prior to reionisation - is independent of den- 
sity. This flattens the temperature-density relation and low- 
ers the value of 7. However, in the SI model, where the 
He II is already in ionisation equilibrium with the UV back- 
ground, the He 11 p hoto- heating rate is instead proportional 
to the IGM density l|Theu ns 2005). The He 11 fraction is pro- 
portional to density (due to higher recombination rates in 
denser regions) and the extra energy injected at z — 3.4 in- 
stead increases 7. As we shall see, these differences will also 
play a small role (relative to To) in the simulated T^g evo- 
lutio n. These values are compared to measurements of 7 
from lSchave et al.l l|2000l ). Although the error bars are large, 
somewhat lower values of 7 are preferred at z < 3.5. 



4.3 Time integration in GADGET-2 and Hydra 

We run the fifth and final model, S3, to quantify the ef- 
fect of the GADGET-2 time integration scheme on our re- 
sults. The study of T02 used hydrodynamical simulations 
perfor med with a modified ve rsion of the P^M - SPH code 
Hydra (|Couchman et al.lll995l ) . Hydra employs a single step 
time integration scheme where all particle positions are ad- 
vanced on the minimum timestep required throughout the 
simulation volume. In contrast, GADGET-2 uses individual 
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timesteps for each particle. This means that particles in high 
density regions, where dynamical timescales are short, have 
timesteps which are orders of magnitude sm aller than part i- 
cle timesteps in the lowest density regions l|Springelll20"05l ). 
This enables efficient use of computational resources and 
substantially improves code performance. However, if a sud- 
den boost to the IGM temperature occurs in the middle 
of a particle timestep, the energy injected will be smeared 
over the timescale corresponding the particle timestep. In 
regions where the dynamical timescale is long, such as the 
low density IGM which dominates the transmission in the 
Lya forest at z = 3, this can delay the impact photo-heating 
on the gas and could potentially affect our numerical results. 

Unfortunately, forcing a single minimum timestep in 
GADGET-2 to test this possibility would be prohibitively 
expensive. We instead impose a stricter upper limit on the 
maximum timestep size for all particles in the S3 model. In 
our first four simulations, the maximum allowed timestep 
for all particles is Aln[l/(1 + z)] = 0.05, although in prac- 
tice this value will vary and can be up to a factor of two 
smaller at any given redshift. This is equivalent to specifying 
the maximum allowed timestep as a fraction of the current 
Hubble time in the simulations. In the S3 model we instead 
impose a maximum timestep of Aln[l/(1 -I- z)] = 0.001 at 
z < 4.2. This choice should provide a useful test of the im- 
pact of timestepping on the gas hydrodynamics during a 
rapid change in the IGM temperature (V. Springel, private 
communication). The thermal and ionisation history of the 
S3 model corresponds to the short-dashed curves in Fig. [S] 
Note the differences between the thermal and ionisation his- 
tories of the SI and S3 models are minimal, and that the 
curves are almost indistinguishable in all three panels. 
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Figure 3. Comparison of Teff measured from our synthetic 
Lya forest spectra to the observational data of FGOSb (open 
diamonds). The error bars correspond to the statistical uncer- 
tainties only. The synthetic data correspond to the models with 
a sharp temperature boost at 2; = 3.4 (SI, solid curve), an 8 per 
cent increase in the free electron fraction in addition to the sharp 
temperature boost (S2, dotted curve) and a stricter upper limit 
on the particle timesteps (S3, dashed curve). The latter model is 
indistinguishable from the SI data. A fourth curve (dot-dashed) 
corresponds to the effective optical depth measured from the SI 
model ignoring the effect of peculiar velocities on the Lya forest. 
As noted by T02, neglecting peculiar velocities raises the absolute 
value of Tcfi by deblending absorption lines. 



5 THE Lya EFFECTIVE OPTICAL DEPTH 

FROM HYDRODYNAMICAL SIMULATIONS 

5.1 Construction of synthetic Lya spectra 

We now turn to analysing the Lya forest effective opti- 
cal depth in each of our simulations. Synthetic Lya forest 
spectra are constructed from each simulation by extract- 
ing randomly selected sight-lines parallel to the simulation 
box boundaries at redshift intervals of Az — 0.05 over the 
range 2 < z < 4. At each redshift a total of 1024 sight-lines 
are extracted, each with 1024 pixels. Every pixel in each 
sight-line has a neutral hydrogen number density nm, tem- 
perature T, peculiar velocity Wpec and Hubble velocity hh 
assoc iated with it. A stan dard SPH interpolation procedure 
{e.g. iTheuns et all 1 19981 ) is used to extract the first three 
of these quantities from the simulation data. In each line of 
sight with A'' pixels, the Lya optical depth in pixel i is then 
given by 

N 

Here bm = {2kBT/muy^'^ is the Doppler parameter, aa = 
4.48 X 10~^* cm^ is the Lya cross-section, SR is the pixel 
width and H{a,x) is the Voigt-Hjerting function (|Hiertinel 
1 19381 ) 



H{a,x) = - dy, (6) 

where x = \vv,{i) - "(j) = "hO') + Wpcc(j), 

a = AaAa/47r6Hi(i), Aa = 6.265 x 10* s~^ is the damp- 
ing constant and Aq = 1215.67 A is the wavelength of 
the Lya transition. We use the analytic a pproximation for 
H{a,x) provided bv lTepper Garcfal (|2006l ). 



5.2 Comparison to the TbS evolution observed by 
FGOSb 

We shall compare our simul ation data to the recent direct 
measurement presented by IFaucher-Giguere et al.l l|2008bl ) 
(hereafter FG08b), who found a feature in the Lya forest 
opacity e volu ti on consistent with the one detected statisti- 
cally by Bernardi et al. (2003). FG08b measured t^s over 
the redshift range 2 < z < 4.2 using a sample of 86 high- 
resolution, high signal-to-noise quasar spectra obtained us- 
ing three different instruments; the ESI and HIRES spectro- 
graphs on Keck and the MIKE spectrograph on Magellan. 

We firstly renormalise the optical depths of the syn- 
thetic spectra by the same constant at every redshift, 
A = 0.87, to approximately match the normalisation of the 
Lya effective optical depth measured by FG08b. For a set 
of spectra with N pixels 
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N 

e--- = (F) = l^e-^-S (7) 

where n is the optical depth in each pixel of the synthetic 
spectra and {F) is the mean observed flux. From Eq. ((2]), this 
renormahsation corresponds to a straightforward rescaUng 
of r_i2 upwards by around 15 per cent if aU other param- 
eters in the simulation remain fixed. We have verified this 
renormahsation has no impact on the shape of the Toff evo- 
lution. 

Fig. [3] shows the comparison between the Tcff evolution 
measured from our synthetic spectra and the observational 
data of FG08b. Note we have yet not processed the synthetic 
spectra to resemble the FG08b data (aside from the over- 
all renormahsation); the synthetic spectra are noiseless and 
have a substantially longer path length than the observa- 
tional data, but nevertheless provide a useful illustration of 
the Tofi evolution in the simulations. The solid, dotted and 
dashed curves display the SI, S2 and S3 models, respectively. 
These data all display the same generic features observed 
in the semi-analytic model; a downturn in t^s due to the 
reheating of the IGM, followed by a slow recovery driven 
mainly by adiabatic cooling. In all instances it is clear the 
simulated data do not exhibit the narrow feature observed 
in the FG08b data. 

Firstly, note the S3 and SI models are again almost in- 
distinguishable, giving us confidence that the adaptive time 
integration scheme used in GADGET-2 will not significantly 
alter our results. Secondly, a small difference between the 
Teff evolution observed in the SI and S2 models is apparent 
at z < 3.4 - recall that the S2 simulation includes the 8 per 
cent increase in the free electron abundance due to sudden 
He II reionisation. As noted in §2, the extra electrons will 
reduce the If ii recombination timescale and thus increase 
the LyQ opacity. In a highly ionised IGM the equilibration 
timescale is icq — frntrac, where /hi = nm/nn is the IGM 
H I fraction and the recombination timescale is 

(l + XHc)A UxIO^kJ V 4 J ■ 

Therefore, for /hi — 10~^'^, the equilibration timescale is 
already very short at mean density, toq ~ 3 x 10* yrs (c/. 
2.9 X 10* yrs between z = 3.4 and z = 3), and the increase 
in /hi (and hence Tcb) following the additional 8 per cent in- 
crease in the electron number density should be very rapid. 
However, we instead find Toff in the S2 model is 3-4 per cent 
lower than the SI model at 3 < 2 < 3.4. In this instance 
the effect of the extra electrons on the Lya opacity is coun- 
tered by the associated flattening of the temperature-density 
relation following sudden, homogeneous He 11 reionisation. 
Recall that most of the Lya forest at z > 3 is dominated 
by transmission from underdense regions in the IGM, which 
become hotter in the S2 model as the temperature-density 
relation flattens. Consequently, including a sudden change in 
the ionisation state of He 11 in addition to a large tempera- 
ture boost increases rather than decreases the magnitude of 
the dip in Tes at z < 3.4. The extra electrons therefore do 
not help in reproducing the Tea data in our model. 

The fourth, dot-dashed curve corresponds to spectra ex- 
tracted from the SI simulation, but this time excluding the 
effect of redshift space distortions induced by peculiar veloc- 
ities in the IGM. As noted by T02, excluding the peculiar 



Table 2. Instrumental resolution (FWHM), pixel sizes and 
signal-to-noise properties adopted for the synthetic Lyo forest 
spectra. The values are based on those reported by FGOSb. 



Instrument 


FWHM 


Pixel size 


S/N per pixel 


HIRES 


6 km 


2 km s-l 


30 


ESI 


33 km s~^ 


11 km s-i 


30 


MIKE 


11 km s^^ 


2 km s-l 


20 



velocities tends to deblend Lya absorption lines, increasing 
Toff. It is also apparent the peculiar velocity field flattens 
the Tefi evolution somewhat; the gradient of the dot-dashed 
curve is steeper in comparison to the SI data including pecu- 
liar velocities, especially at z < 3. These differences decrease 
towards lower redshift, indicating line blending becomes less 
widespread as the line number density and Lya opacity both 
fall. However, in contrast to T02 (see their fig. 2), we find 
these redshift space distortions are insufficient to produce 
the rapid upturn in by z = 2.9 seen in the observa- 
tional data. The redshift space distortions induced by pecu- 
liar velocities provide no change in the generic shape of the 
Toff evolution with redshift. We have verified that this also 
holds for the S3 model with the stricter timestep limit. This 
result will be considered in more detail in !j6l 

However, to make a fair comparison to the FG08b data, 
our synthetic spectra must be processed to resemble their 
observational data set as closely as possible. Our renor- 
malised synthetic spectra are therefore also convolved with 
a Gaussian instrument profile and resampled onto pixels of 
the required size. Gaussian distributed random noise is then 
added to each pixel. The parameters used for this proce- 
dure are summarised in Table [2] corresponding to the three 
instruments used in the FG08b data set. A random sub- 
sample of synthetic HIRES, ESI and MIKE spectra with a 
total path length corresponding to the values displayed in 
fig. 2 of FGOSb are then drawn from our simulated data set 
in intervals of Az = 0.2 between 2 < z < 4. 

The results are displayed as the open diamonds in Fig.|4] 
for the SI (left panel), El (central panel) and Nl (right 
panel) models. Following FGOSb, the error bars (statistical 
only) correspond to the standard error of the mean, com- 
puted by subdividing the synthetic data into chunks 3 proper 
Mpc in length. The solid curves in each panel correspond to 
the underlying Tcff evolution measured from the unprocessed 
synthetic Lya spectra. Our results are compared to the best 
fitting curve to the FGOSb data, displayed by the dashed 
line in each panel. 

As expected, the synthetic data in Fig. |4] now exhibit 
some additional scatter due to the variation in cosmic struc- 
ture probed by the random sight-lines, but in all instances 
are within la of the underlying Tcff evolution. The SI data 
in the left-hand panel again exhibit the generic behaviour 
seen in the semi-analytical model, and are inconsistent with 
the narrow feature observed by FGOSb. For comparison, the 
central panel displays the data obtained fro m the El model 
with an extended period of heating {e.g. iMcO.uinn et al.l 
l2009h . while the right hand panel displays the Nl model 
with no additional heating from He 11 reionisation. In both 
instances t^b smoothly evolves with redshift. This result is 
in agreement with the predictions for Tcfr from the radiative 



© 0000 RAS, MNRAS 000, 000-000 



Evolution of the Lya optical depth 9 




0.1 1 1 

2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 



z z z 

Figure 4. A comparison of the TcS evolution measured from synthetic Lya spectra drawn from three of our hydrodynamical simulations 
to the best fitting curve from FGOSb (dashed curves). The solid curves in each panel display r^ff measured from the unprocessed synthetic 
Lya spectra, while the open diamonds with error bars correspond to r^ff measured in bins of width Az = 0.2 after the spectra are processed 
to closely resemble the FGOSb data set. The error bars correspond to the standard error of the mean. Left: The SI model. Centre: The El 
model. Right: The Nl model. 



transfer simulations of iMcQuinn eT&il l|2009l '): an extended 
period of reheating during He ii reionisation will induce no 
sharp features in the t^s evolution due to changes in the 
IGM temperature alone. 

Our detailed simulations therefore indicate that hydro- 
dynamical effects in the IGM following rapid reheating (but 
see rBolton et al.ll2009l : iMcQuinn et al.ll2009h wiU not aid in 
reproducing the sharp dip in the Lya opacity observed by 
FGOSb. Instead, our results suggest any recovery in TcS fol- 
lowing sudden reheating will be driven primarily by adia- 
batic cooling in the low density IGM, consistent with the 
behaviour predicted in our semi-analytical model. 



6 HYDRODYNAMICS AND REDSHIFT 
SPACE DISTORTIONS FOLLOWING 
SUDDEN REHEATING 

In the last section we demonstrated our simulations were 
unable to reproduce the narrow feature observed by FGOSb. 
We now consider this result in more detail by examining the 
impact of hydrodynamical effects and redshift space distor- 
tions on the Lya effective optical depth evolution. 



6.1 The gas density distribution 

We begin by briefly discussing the effect of sudden reheat- 
ing on the IGM density distribution. Following an injec- 
tion of energy into the IGM during reionisation, the sub- 
sequent increase in gas pressure will smooth the IGM den- 
sity distribution over s cales correspondi ng to the local Jeans 
len gth jSc hayc 2001; P awlik et al.[|2009l '). However, as noted 
bv lGnedin fc Hui (199§ ), this Jeans smoothing will not oc- 
cur instantaneously. The low density gas responsible for the 
Lya forest absorption will expand on a sound- crossing time, 
which can be considerable - comparable to the Hubble time 
- in low density systems. Indeed, hydrostatic equilibrium 
is only restored once tsc ~ tdyn, where isc and t^yn are 
the sound cr ossing and dynamical timescales respectively 
(|Schavell200il ). The dynamical time fdyn in the low density 



IGM is around the Hubble time, as is evident from the Fried- 
mann equation, = SirGp. 

This can be appreciated in the comparison between 
sight-line data drawn from the SI (solid curves) and Nl (dot- 
ted curves) simulations displayed in Fig. [S] The left panel 
displays a subsection of a single sight-line through the simu- 
lations a,t z — 3, while the right panel shows the same sight- 
line ai z — 2. From top to bottom, the normalised gas den- 
sities, the H I fractions, the gas temperatures, the peculiar 
velocities and the resulting Lya forest spectra, computed in 
three different ways, are displayed. The spectra in the third 
row from bottom are computed using Eq. ((5} and include 
the effect of redshift space distortions due to the peculiar 
velocity field. In the next row down, however, the Nl spec- 
trum has been recomputed by fixing the gas temperature in 
the Nl model to be equal to the SI values and then rescal- 
ing the Nl H i fraction, such that n^j = nm(?si/TNi)~'''^. 
This removes most of differences in the spectra attributable 
to Doppler broadening and the H i fraction, although note 
the r~*^ '^ scaling for the recombination coefficient is not 
exact and small differences will remain Finally, in the bot- 
tom row the SI and temperature adjusted Nl spectra both 
exclude the effect of the peculiar velocity field. 

It is clear that differences in the Lya spectra due to 
the gas density alone are very small soon after reheating. 
The density distribution at z; = 3 has only been slightly 
smoothed around the peak of the overdensity located at 
6.9h~^ Mpc. In contrast, pressure smoothing of the gas dis- 
tribution is clearly apparent hy z — 2, as there has been 
sufficient time for the gas distribution to be altered by the 
increased pressure in the high density regions. It is precisely 
because of the finite timescale required to achieve hydro- 
static equilibrium that changes in the gas density distribu- 
tion due to Jeans smoothing will have a negligible impact 



* In the left hand panel of Fig.[5]prior to this rescaling, the larger 
difference between the SI and Nl line profiles at 5.Sh~^ Mpc on 
comparison to the broader lines at 6.7h~^ Mpc is because these 
absorption features lie on the linear and logarithmic parts of the 
curve of growth, respectively. 
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Figure 5. Comparison between various IGM properties along a subsection of a sight-line extracted from identical locations in the SI (solid 
curves) and Nl simulations (dotted curves) at 2 = 3 (left panel) and 2 = 2 (right panel). From top to bottom, the normalised gas densities, 
the H I fractions, the gas temperatures, the peculiar velocities and the resulting Lya forest spectra, computed in three different ways, are 
displayed. The spectra in the third row from bottom are computed using Eq. JSj and include the effect of the peculiar velocity field. In 
the next row, however, the Nl spectrum has been recomputed by fixing the gas temperature in the Nl model to be equal to the SI values 
and then rescaling the Nl H I fraction, such that nj^j = nHi(T'si/TNi)~'''^. This removes most of differences in the spectra attributable 
to Doppler broadening and the H I fraction. The SI and temperature adjusted Nl spectra both exclude the effect of the peculiar velocity 
field in the bottom row. 



on the Teff evolution by z = 3 if the IGM is suddenly heated 
aX z = 3.4. 

6.2 Peculiar velocity gradients 

However, the finite time required to change the gas distri- 
bution does not exclude rapid changes in the peculiar veloc- 
ity gradients responsible for the Jeans smoothing. Redshift 
space distortions associated with these g radients may then 
induce changes in the Lya line profiles (|Brvan et al.lll9S)3 : 
euns et al.ll2000l) . T02 appealed to these peculiar velocity 
gr adients to explai n the r ecovery in the Toff feature observed 
by iBernardi et al] (|2003l ). On examination of the peculiar 
velocity field in the left panel of Fig. [S] the SI model does 



indeed exhibit a positive peculiar velocity gradient in the 
centre of the overdensity at Q.^hT^ Mpc, as it must if the 
gas is to expand. The increased gas pressure has reversed the 
collapse of the overdensity which is occurring in the colder 
Nl model. 

This effect of gas temperature on the peculiar veloc- 
ity gradients is displayed in more detail in the left hand 
panel Fig. (6] where the probability distribution of the pe- 
culiar velocity derivative with respect to the Hubble veloc- 
ity, in all 1024 synthetic sight-lines is shown for the 
SI (solid curve), El (dotted curve) and Nl (dashed curve) 
simulations. The distribution is shown for overden se gas 
with A > 3 only (c/. fig. 3 in iTheuns et al]l2000l ). The 
shift in the probability distribution towards more positive 
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Figure 6. Left: Probability distribution of the peculiar velocity derivative with respect to the Hubble velocity, tjjj, in the SI (solid curve), 
El (dotted curve) and Nl (dashed curve) simulations at z = 3. The probability distribution is for gas with A > 3 only. The shift in the 
probability distribution towards more positive gradients for progressively hotter models is associated with expansion in the higher density 
regions in the simulations. Right: The probability distribution function of the difference between the peculiar velocity fields in the SI and 
Nl models in all regions where there is a positive velocity gradient in the SI model {i.e regions which are expanding). The solid and dotted 
curves display the distributions for all gas with A > 3 and A > 10, respectively. The largest values of Aiipoc correspond to the densest 
regions in the simulations. 



gradients for progressively hotter models is due to expan- 
sion in the hotter, high density regions which have become 
overpressurised wi t h respect to their cooler surroundings 
iBrvan et aLlll999l : iTheuns et aLll2000l '). The fact that the 
probability distributions look qualitatively s imilar to the 
data presented in fig. 3 of iTheuns et al.1 (|2000l ) IS encourag- 
ing, and indicates the higher temperatures in our GADGET- 
2 sim ulations have a simil ar impact on the peculiar velocity 
field. ITheuns et all (|2000l ) used simulations performed with 
Hydra in a 2.5/i~^ comoving Mpc box with a gas particle 
masses of 1.14 X 10*^ Mq (similar to our mass resolution). 
Note the larger velocity gradients present in our simula- 
tions are most likely due to the additional large scale power 
present in our significantly bigger simulation boxes {lbh~^ 
comoving Mpc). 

In the example sight-line in Fig. O the maximum differ- 
ence between the peculiar velocities in the SI and Nl models 
is ~ 10 km s~^ in the centre of the overdensity. Although 
this difference indeed has some impact on the line profile, 
the change in the peculiar velocity field is small on compar- 
ison to the line width, and it is not enough to significantly 
alter the broad Lya line shown in the lower panels. Once 
differences attributable to the different temperatures of the 
simulations have been scaled out of the spectra, the SI and 
Nl models produce very similar absorption line profiles. 

This anecdotal evidence is displayed more quantita- 
tively in the right hand panel of Fig. [6] where we plot the 
probability distribution of the difference between the pecu- 
liar velocities in the SI and Nl simulations for all pixels 
with positive peculiar velocity gradients and A > 3 (solid 
curve) or A > 10 (dotted curve) in the SI data. In the ma- 
jority of these regions there is in fact only a small change 
in the peculiar velocity field, with the largest differences 



(~ 10 km s~^) associated with the rarer, high density peaks 
which produce broad lines like the one shown in Fig. [5] Since 
an absolute change in the peculiar velocity of a few km s~^ 
is s mall in com parison to typical line widths of 20 km s~^ 
(e.Q. lKim et al. 2002), this explains why the impact on t^s is 
correspondingly small. Instead, it is the instantaneous tem- 
perature of the IGM which primarily influences the TeS evo- 
lution in our simulations. 



These flndings differ from those of T02, who found 
changes in the peculiar velocity gradients were partially 
responsib le for inducing the recovery in Tch observed by 
iBernardi et al.l (|2003l ). The discrepancy between these re- 
sults may be due to differences between the numerical meth- 
ods used. We have tested the impact of timestepping on sim- 
ulations, and we find our results are robust in this respect. 
However, we cannot be absolutely certain that other numer- 
ical effects do not play a role, and ultimately we can only 
speculate on the origin of this difference. Ideally, an inde- 
pendent numerical study is required to verify or refute our 
claims. In agreement with T02, however, we do indeed find 
that line blending due to the peculiar velocity field lowers 
Toff, and that higher temperatures steepen the peculiar ve- 
locity gradients in overdense regions, giving us confidence 
that our simulations are at least broadly consistent with 
T02. We find these effects are nevertheless insufficient to 
reproduce the narrow feature observed in the Tbs evolution, 
suggesting that a sudden increase in the IGM temperature at 
z = 3.4 following He ii reionisation is unable to adequately 
explain the Tch data. 
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r ON THE POSSIBLE ORIGIN OF THE 
OBSERVED Tcff FEATURE 

Thus far we have estabhshed that even in the presence of a 
sudden increase in the IGM temperature at z = 3.4 (but 
see iBohon ct al. 2009; McQuinn ct al. 20Qii), the narrow 
dip observed in the Lya effective optical depth by FGOSb 
cannot be reproduced in our hydrodynamical simulations 
of the Lya forest. However, before proceeding it is worth 
briefly emphasising that the purported "narrowness" of the 
Teff feature is a somewhat model-dependent statement. For 
instance, if we consider the Tcff — z plane, FGOSb charac- 
terised the width of the feature by fitting a power-law and 
a Gaussian "bump" . This significance of the bump then de- 
pends on how good an approximation a power law is to the 
underlying evolution of the IGM opacity excluding the ef- 
fects associated with He ii reionisation. In our simulations 
without sudden heating, t^s does indeed evolve smoothly 
with redshift, indicating this should be a reasonable approx- 
imation if the redshift evolution of the ionising background 
is gradual. In contrast, the FGOSb feature stands out visu- 
ally, and more objective measures we have applied to the 
data {e.g., a regularised derivative of the data as used in 
edge- finding algorithms) tend to confirm that impression. 
The fact that it is seen in at least three independent data 
sets further hints that it is not a data artifact. 

Additionally, a decomposition into the required evolu- 
tion in the ionising background (as in the next section) will 
also depend on the assumed IGM thermal evolution. The 
Lya opacity scales as t oc r~'^'^/rHi; the effect of gas tem- 
perature and the photoionisation rate on TeS are impossi- 
ble to disentangle without independent estimates for one 
or the other. In this work we have demonstrated that sud- 
den heating on its own cannot reproduce the Tcff feature. 
However, a sudden global heating event followed by a grad- 
ual downturn in the photo-ionisation rate at 2; < 3.2 could 
still reproduce the data. On the other hand, there is good 
reason to suppose that such a sudden, global boost to the 
IGM temperature is unlikely in the first place. Our pre- 
vious work has demonstrated that He 11 reionisation likely 
produces too little heating to produce a subs tantial opac- 
ity d ecrease over the short timescale required (|Bolton et al.l 
I2OO9I ). In the next section we therefore assume the tempera- 
ture boost during He 11 reionisation occurs over an extended 
period of time (our El model). This model also resembles the 
res ults from recent radiat ive transfer simulations performed 
bv lMcQ uinn et al] l|2009l ). We therefore dispense with rapid 
changes in the temperature entirely, and now turn to dis- 
cuss the remaining possibility for the origin of the Toff dip: a 
narrow peak in the metagalactic hydrogen photo-ionisation 
rate. 



7.1 The hydrogen photo-ionisation rate required 
by the FGOSb Tcff evolution 

Many studies have used simulations of the Lya for- 
est, combined with measurements of r^ff , to place con- 
straints on the metagalactic hydrog en photo-ion i sation 
rate, F- 



Fhi/IO"^^ s~^ (e.g. iRauch et al.1 



199 



200 



McDonald fc Miralda-Escudell200ll: I Cen fc McDonald! 
Schavc ct al. 2003; Meiksin & Wh iti |2004| : 
20041 : iBolton et al.ll2005l : iJena et al.ll2005l '). We use the same 
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Figure 7. The solid curve corresponds to the hydrogen photo- 
ionisation rate, r_i2 = Fhi/IO"^'^ required to reproduce 
the fit to Tcff evolution presented by FGOSb. The narrow peak in 
r_i2 at 2 = 3.2 corresponds to the dip in the FGOSb Teff evo- 
lution. The dashed curve displays the ex pected contribu t ion o f 
quasars to the total F_i2, based on the iHopkins et al. 
model for the quasar luminosity function, and the dotted line 
corresponds to the difference between the total F_i2 required by 
Tcff and the quasar only contribution. To give an indication of 
the statistical uncertainty in this peak, the open diamonds with 
error bars are derived from the FGOSb Teff data points including 
their la statistical errors. 



procedure in this work to estimate r_i2 from our El sim- 
ulation, which was constructed to have a similar thermal 
history at mean density to t he rec ent radiative transfer sim- 
ulations of iMcQuinn et al.l (|2009l). We use Eq. ^ and the 
best filQ to Teff obtained by FGOSb for this purpose. Fol- 
lowing from Eq. ([2]), the photo-ionisation rate which repro- 
duces the Teff fit is given by r_i2(z) = rHM{z)/A{z), where 
Lhm{z) is the photo-ionisation rate from the HMOl UVB 
model divided by 10~^^ s~^. 

The result of this procedure is displayed as the solid 
curve in Fig.[7l A gradual increase in r-_i2 towards lower red- 
shift is required, and the prominent bump centred ai z — 3.2 
is responsible for the narrow Tcff feature. The open diamonds 
with error bars are derived from the FGOSb Tcff data points 
including their la statistical errors; these give an indication 
of the statistical uncertainty in this peak. 

The dashed curve shows the expected contribution 
from quasars to r-12 . We c ompute this using the re- 
cent 



Madau fc Haardd (l2009h parameterisation for the 



iHopkins et all 1120071 1 comoving quasar emissivity at the 



H I Lyman limit, ui,, ass uming a power law spectrum with 
= eL(j^/i^)"^" (e.g. iTelfer et all l2002[ l. The expected 
contribution to the photo-ionisation rate from quasars is 
then 



® Specifically, we use the best fit FGOSb present for their Teff mea- 
surements in redshift bins of width Az = 0.2 when using the 



ISchave et al. I ll2003h metal correction. 
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(1 + ^) 

10-12 s 



(9) 



6.3 X 



where the photo-ionisation cross-section 
10~^* cm'^ (v and the mean free path = 
Xh{i'/i/LY'^; the latter assumes the H i column densi t y dis- 
tribution f(NTii,z) PC N-^-^ {e.g. iPetitiean et~ai1 Il993l : 
lMiralda-Escudil2003l ). We take Al = 50[{1 + z)/4]-^ proper 
Mpc, which is towards the lower end ofthe range recently ad- 
vocated bv iFaucher-Giguere et ahl (|2008al ). Note, however, 
that the expected contribution to the photo-ionisation rate 
from quasars depends on the uncertain quasar luminosity 
function (particularly the faint end slope), mean free path 
and ionising spectrum, as well as radiative transfer effects 
(for in stance, reprocessed rad iation can contribute signifi- 
cantly; iHaardt fc Madaull 19961 ) . 

The dotted curve displays the difference between the to- 
tal photo-ionisation rate and the contribution from quasars 
alone, r'_i2 = r_i2 — r5i2- This curve remains almost flat 
for 2 < z < 4, with r'_i2 — 0.55, aside from the narrow 
bump {Az ~ 0.4) which peaks with r'_i2 — 0.75 at z — 3.2. 
Thus, in order to reproduce the observed Tcb, we require a 
boost of around 35 — 40 per cent in r'_i2 over a short redshift 
interval. 

Finally, note that although the FGOSb statistical er- 
rors are relatively small, the normalisation of r_i2 is still 
somewhat uncertain. In particular, although the r_i2 values 
we derive from our El simulation are consiste nt with other 
estimates from hydrody n amical simulation s (jTvtler et al.l 
l2004l : iBolton et ahl l2005l : Ijena et all |2005| ). they are sys- 
tematically higher by up to a facto r of two in comparison 
to the recent estimates presented bv IFaucher-Giguere et ahl 
l^2008a^ using the sarn e Teg da ta. This is partially because 
IFaucher-Giguere et ahl |20083) do not model the Lya for- 
est in detail, instead obtaining analytical estimates for r_i2 
using the FGP A combined with the I GM d ensity distribu- 
tion derived bv lMiralda-Escude et al.l (|200G| ') (see Eqs.[2]and 
[3|. Systematic uncertainties on r_i2 due to other parame- 
ters, such as the IGM t emperature which i s still poorly con- 
strained, are also large (|Bolton et aLlfioOsi ). Thus, although 
the overall normalisation remains somewhat uncertain, the 
shape of these curves should be fairly robust, and the re- 
quirement for a peak in r_i2 remains so lon g as the volume 
weigh ted IGM temperature evolves slowly ( McQuinn et al.l 
I2009D . 



7.2 Possible causes of a peak in the hydrogen 
photo-ionisation rate 

We have found that, even if the IGM temperature increases 
suddenly, we cannot reproduce a narrow feature in the opti- 
cal depthQ Thus a narrow peak at 2; = 3.2 in the otherwise 
approximately constant value of r'_i2 = r_i2 — r5i2 from 
2 < z < 4 is instead required to reproduce the t^b evolution 
measured by FGOSb. This could be associated with a sharp 
modulation in the ionising emissivity or the mean free path 



° The exception to this is if the photo-ionisation rate subse- 
quently drops at 2 < 3.2 following a sudden temperature in- 
crease. However, as discussed earlier, recent studies indicate that 
a rapid global temperature boost is un likely in the first place 
llBolton et al.ll2009l: iMcQuinn et al.ll2009l) . 



of ionising photons, or indeed a combination of both effects 
(see eq. [9)l . However, while appealing to the stellar contri- 
bution to the photo-ionisation rate for an increase in the 
emissivity is in principle acceptable given the uncertainties 
involved in derivi ng the stellar ionising rate from observa- 
tions at z > 3 (e.q.lMadau et al.lll999l : lMiralda-Esc"udil2003l : 



iBohon fc Haehneltll2007bl v" there is no apparent reason why 
such a narrow peak should occur at z = 3.2. 

Instead, we suggest it is more likely that any modu- 
lation of the hydrogen photo-ionisation rate is associated 
with dense He 11 and H i Lyman limit systems. Unlike the 
low density IGM, such systems have dynamical, cooling and 
recombination times which are comparable or shorter than 
the ~ 10* years associated with the width of the TcB fea- 
ture. They can modulate the hydrogen ionising background, 
either by changing the hydrogen ionising emissivity (by the 
reprocessing of He 11 ionising photons into He 11 Lyman se- 
ries, Balmer, or two-photon emission, all of which can ionise 
hydrogen), or the mean free path of ionising photons, since 
their opacity will be altered by changes in their size or tem- 
peraturcQ- Furthermore, a sharp change in the emissivity or 
opacity will then be imprinted on the ionising background 
on a timescale comparable to the mean free time of an ion- 
ising photon, or Az ~ {dN/dz)~^[r{l3 - 1)]"^ ~ 0.24 where 
dN/dz = 3.3[ (1 -|- z)/5\^'^ is the abundan ce of H i Lyman 
limit systems (I Storrie-Lombardi et al ]ll994'), and we assume 
dN/dNui oc N^f with /? = 1.5. This is certainly well within 
the range required to explain the feature. We now proceed 
to give some brief examples of these effects. 

• Helium recombination radiation 

A large fraction of the hydrogen ionising background may 
be radiation from higher freque ncies that is rep r ocess ed by 
dense systems. For instance, iHaardt fc Madaul (|l99(3l ) find 
that ~ 40 per cent of Phi comes from re processed radia- 
tion at 2 ~ 3; similarly (|Fardal et al.|[l99^ ) find that about 
~ 20 per cent of Fhi comes from reprocessed radiation. 
This has two components: recombinations of hydrogen to 
the ground state, and reprocessing of He 11 ionising photons 
into He 11 Lyman series, Balmer, or two-photon emission, 
all of which can ionise hydrogen. The latter obviously under- 
goes rapid evolution during the process of He 11 reionisation, 
particularly toward the tail end of reionisation when the 
mean ionising background can rise rapidly (although there 
could be large fluctu ations in the backg round throughout the 
reionisation process; iFurlanettdboOsi ). On the other hand, 
recent calculations indicate the contribution of recombina- 
tion radiation to th e UV background could be s maller than 
previous estimates (|Faucher-Giguere et al.|[2009l ). 

Another potentially important effect comes from the 
reprocessin g of He 11 Lym an series photons in an optically 
thick IGM (|Madau fc Ha ardt 2009). Lyman series photons 
between Ly/3 (at E — 3.56 Ry) and the Lyman limit 
(4 Ry) are degraded to He 11 Lya photons, Balmer or 
lower frequency radiation. The magnitude of this effect 
depends strongly on the opacity of the IGM, which of 
course evolves rapidly during He 11 reionisation. While this 
may not appear to be particularly significant (since such 
photons can already ionise hydrogen), the degradation to 



^ The heating and expansion of Lyman limit systems was also 



i ne neatmg ana expansion or ijyman r 
briefly discussed bv lMcQuinn et al.l ||2009| ') 
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lower energies implies that the hydrogen photo-ionisation 
rate can be significantly modulated. In particular, selection 
rules forbid Ly/3 photons from being converted into Lya 
photons, and the reduction of two-photon emission from 
Ly/3 reprocessing at the end of He ii reionisation could 
reduce the hydrogen ionisation rate, helping to produce the 
downturn in the photo-ionisation rate. 

• The size of Lyman limit systems 

In this work, we have thus far argued that hydrody- 
namic effects associated with the heating of the IGM during 
He II reionisation have little effect on TeS- While this is true 
of the low-density gas responsible for most of the absorption 
in the Lya forest, it is not true of higher density systems 
which have smaller dynamical timescales, fdyn oc A"^''^. 
Since H i Lyman limit systems (LLSs) dominate the opacity 
of the post-reionisation IGM, a sudden change in the sizes 
of these systems can then modulate the mean free path of 
H I ionising photons. 

We may consider this as follows. At the end of 
He II reionisation, although most of the IGM has been re- 
heated, self-shielded He ii LLSs, in which the helium is still 
predominantly in the form of He ii, will still remain. Self- 
shielding to He II ionising photons occurs for He ii column 
densities A^hcII > 6.7 x 10^ cm~ . Assuming the size of 
an a bsorber with overdensity A is the local Jeans length 
jSchave 2001'). this column density corresponds to a charac- 
teristic overdensity 

where r_i4 — FhoIi/IO^^^s^^ is the He ii photo-ionisation 
rate. The H i column d ensity of an absorber with overden- 
sity A is l|Schavdl200il ) 

a3/2 / \ -1/5 /I I _s 9/2 

iV„,.3.8xI0-W^(^) (1±£) .(U) 

Self-shielding to He ii ionising photons therefore hap- 
pens at an Hi column density of A'hi ~ 1.6 x 
10^^ cm~^(r_i4/r_i2), well short of the column density of 
A^'hi ~ 1.6 X lO^'^cm"^ associated with H i LLSs. Thus, at 
face value it may seem that He ii reionisation cannot sub- 
stantially affect H I LLSs, and hence the mean free path of 
H I ionising photons. However, He ii LLSs are also likely to 
be associated with the less dense outer regions of H i LLSs. 
A strong increase in external pressure due to the reheat- 
ing of the IGM during He ii reionisation then means that 
the system will no longer be in hydrostatic equilibrium. The 
pressure gradient compresses the gaslf] until pressure equi- 
librium is once again restored (a similar me chanism ha s been 
invoked for globular cluster formation, e.g. Gen 2Q0li) . As a 
result, H i LLSs could decrease in size, and increase in den- 
sity. Alternatively, if hard photons can penetrate and hea t 
the gas in the high density regions (e.g. iBolton et aDl2009l l. 
the H I LLS could instead expand. 

For simplicity, let us for the moment assume that any 
change in density occurs isothermally (see below for more 

* Note that in our optically thin hydrodynamical simulations the 
self-shielding of dense systems is not modelled. This effect is there- 
fore absent in our simulation data. 



discussion). If the absorber expands or contracts to a new 
overdensity A — > /A, then (since riHi oc A^ and R oc 
A"^''"^), A^Hi c>c f^^^ , while the cross-section of the absorber 
cr (X f~'^^^ . Assuming that / is independent of A (the opac- 
ity is dominated by H i LLSs, so the result will be most af- 
fected by / for such systems) , the mean free path of ionising 
photons will be altered by a factor 

A ^ !Z--'-'[l-eM-r)dr] 

Al /J"^/-2/3^-i.5 [1 „exp(-/V3^)d^]' ^ 

where the lower limit r ~ 0.1 is the H i optical depth at 
which systems self-shield from He ii ionising radiation. This 
is a non-monotonic function with a minimum at / ~ 0.3; 
however, it does not exceed unity until / > 1. Thus, if 
systems compress, the hydrogen photo-ionisation rate falls; 
if they expand, it increases. However, more quantitative 
exploration of this possibility requires careful simulation 
of the LLS with coupled hydrodynamics and radiative 
transfer, particularly since the equilibrium temperature 
also strongly aflects opacity. In addition, a caveat to this 
argument is that although the dynamical timescale of 
an individual Lyman limit system is comparable to that 
required for the modulation of Fhi, a globally averaged 
change in the mean free path still requires the size of all the 
Lyman limit systems to change over a short time interval. 

• The temperature of Lyman limit systems 
Yet another possibility is that the temperature of 
H I LLSs itself changes during He ii reionisation, altering 
their opacity. A rapid increase in the injection of photons 
at 3 Ry from reprocessed radiation implies that the ionising 
background hardens significantly, which could result in a 
significant change in the equilibrium temperature. This 
could be further modulated by the hydrodynamic effects 
mentioned above. Note that the equilibrium temperature of 
LLSs is a non-trivial function of density in photo-ionisation 
equilibrium, particularly in the presence of metal line cool- 
ing (Wiersma et al. 200t^). Generally, as a system becomes 
denser, the equilibrium temperature falls, since the effi- 
ciency of cooling i ncreases as the i onisa tion parameter falls 
dEfstathioul 1 19921 : Iwiersma et al] 120091 '). This exacerbates 
the increase in opacity due to the density increase. Non- 
monotonic evolution of the temperature of the LLSs (which 
have short cooling times since they cool radiatively, rather 
than adiabatically) could thus modulate their opacity, as 
well as the mean free path and the H i photo-ionisation rate. 



8 CONCLUSIONS 

We have used a semi-analytic model of inhomogeneous 
He II reionisation and high resolution hydrodynamical sim- 
ulations of the Lya forest to investigate the impact of sud- 
den reheating on the evolution of the Lya forest effective 
optical depth. Our semi-analytic model indicates that any 
injection of energy into the IGM during inhomogeneous 
He II reionisation will produce a well understood and generic 
evolution in Tcff, where a reduction in the opacity from 
2: = 4 to 2: = 3 is followed by a gradual, monotonic re- 
covery driven largely by adiabatic cooling in the low den- 
sity IGM. This behaviour is inconsistent with the narrow 
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dip {/S.Z = 0.4) of around 10 per cent in Tefi at 2 = 3.2 
which has now been detected by three independent observa- 
tional st udies f'Bcrnardi et al.l [20031 : iDall'Aglio eralll2008l : 
iFaucher-Giguere ct al...2008bf K 

However, our semi-analytic model does not include a 
detailed reconstruction of the Lya forest. We therefore also 
analyse five high resolution hydrodynamical simulations of 
the IGM to investigate the effect of various thermal histo- 
ries on the Lya effective optical depth. We find that sud- 
den reheating ai z = 3.4 results in a sharp decrease in 
Teff, although we note that such a large, sudden increase in 
the IGM temperature is nevertheless unlikely to occur over 
the e ntire IGM at once (|Bolton et al.ll2009l : iMcQuinn et al.1 
l200g| ). This assumption may nevertheless be appropriate in 
localised regions around quasars with hard spectra, and in 
small volumes similar to the box size of our hydrodynamical 
simulations. 
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Although the assumption of sudden, homogeneous re- 
heating does indeed successfully reproduce the initiation of 
the Teff dip observed by FGOSb, our simulations are still un- 
able to account for the rapid recovery of the narrow dip in 
the Teff evolution hy z = 2.9. The effect of the 8 per cent in- 
crease in the free extra fraction following He 11 reionisation 
on Teff is small, and any resulting increase in Toff is instead 
countered by the simultaneous flattening of the power-law 
temperature-density relation during instantaneous, homo- 
geneous He ft reionisation at z = 3.4. Redshift space dis- 
tortions in the Lya forest attributable to the response of 
the gas to the extra energy injected into the IGM are also 
unable to account for the observed recovery of Tcb- We find 
that sudden reheating does indeed produce larger, more pos- 
itive peculiar velocity gradients corresponding to regions of 
newly expanding gas, but the absolute changes in the ve- 
locity field are generally small in comparison to typical line 
widths, except in the most overdense regions in the simu- 
lation. Such regions are, however, rare and thus contribute 
little to the average opacity. We have tested the robustness 
of this result with respect to the time integration scheme 
employed in GADGET-2, and we find this has little impact 
on our numerical results. In contrast, if He 11 reionisation 
is an ext ended process, t h en in agreement with the recent 
study by iMcQuinn etahl (|2009l ) we find the Lya effective 
optical depth will evolve smoothly with redshift. 

As a consequence, we must instead appeal to a narrow 
peak in the metagalactic hydrogen photo-ionisation rate at 
z — 3.2 to reproduce the Toft feature in our simulations. This 
could potentially be modulated by Lyman limit systems, 
which have recombination, cooling and dynamical times 
comparable to or less than the timescale associated with 
the width of the feature. In particular, we suggest that ra- 
diative transfer effects from He 11 reionisation itself could be 
responsible, either by altering the emissivity of reprocessed 
He ft recombination photons, or by changing the opacity of 
H f Lyman limit systems and hence the mean free path of 
ionising photons. However, further detailed investigation is 
still required to establish the origin of this intriguing feature 
in the redshift evolution of the Lya forest opacity. 
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